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We study cosmological self-reproduction in models of inflation driven by a scalar field (j> with a 
noncanonical kinetic term (k- inflation). We develop a general criterion for the existence of attractors 
and establish conditions selecting a class of fc-inflation models that admit a unique attractor solution. 
We then consider quantum fluctuations on the attractor background. We show that the correlation 
length of the fluctuations is of order CsH" 1 , where c s is the speed of sound. By computing the 
magnitude of field fluctuations, we determine the coefficients of Fokker-Planck equations describing 
the probability distribution of the spatially averaged field (j>. The field fluctuations are generally 
large in the inflationary attractor regime; hence, eternal self-reproduction is a generic feature of k- 
inflation. This is established more formally by demonstrating the existence of stationary solutions of 
the relevant FP equations. We also show that there exists a (model-dependent) range (j}R < <f> < 0max 
within which large fluctuations are likely to drive the field towards the upper boundary <f> = ma x, 
where the semiclassical consideration breaks down. An exit from inflation into reheating without 
reaching ma x will occur almost surely (with probability 1) only if the initial value of <f> is below <f>R- 
In this way, strong self-reproduction effects constrain models of fc-inflation. 



I. INTRODUCTION 

The paradigm of fc-inflation pj assumes that the energy 
density in the early universe is dominated by a scalar field 
<t> with an effective Lagrangian p(X, (j>) consisting solely 
of a (noncanonical) kinetic term, 



p{X,<j>)=K{(j>)L(X), X 



(1) 



where K(<f>) and L(X) are functions determined by 
the underlying fundamental theory. The effective La- 
grangian is assumed to be valid in a sufficiently wide 
range of values of <f> and d^cj), as long as the energy den- 
sity of the field <f> is below a suitable scale (e.g. the Planck 
scale). Unlike the "traditional" models of inflation, where 
the dominant energy density is due to a potential V{4>), 
an accelerated expansion in fc-inflation is driven by the 
kinetic energy of the scalar field. 1 The form (JIJ of the La- 
grangian is sufficiently general to cover many interesting 
cases. For instance, Lagrangians of the form 



p(X, 0) = ^i(0)X ni + K 2 {4,)X r 



(2) 



where n\,rt2 are fixed numbers, are reduced to the La- 
grangian after a suitable redefinition of the field <j>. 

The assumed initial conditions for fc-inflation consist 
of the field 4> having small spatial gradients within a suf- 
ficiently large initial region, and values of (f> and the time 
derivative <fi within suitable (and wide) ranges. These 



1 An effective scalar field <j> with noncanonical kinetic terms is also 
used in models of k-essence l[3Qj and more recently e.g. 010 
where the Lagrangian p(X, <j>) is arranged to track the equation 
of state of other dominant matter components and to become 
dominant only at late times. In the present paper, we do not 
consider such models of fc-essence but assume that the energy 
density of the field <j> is dominant during inflation. 



initial conditions guarantee that the energy density due 
to the field <fi is positive and dominates other forms of en- 
ergy in the universe. The evolving homogeneous field <f>(t) 
rapidly approaches an attractor solution, 4>*{t), which 
drives inflation. The presence of the attractor makes a 
model largely insensitive to initial conditions, since the 
attractor solution is unique and is approached exponen- 
tially quickly. In the present paper we develop a general 
definition of attractor solutions and demonstrate their 
existence in a wide range of models of fc-inflation. 

The back-reaction of quantum fluctuations modifies 
this classical picture of cosmological evolution, adding 
random "jumps" to the classical "drift" along the attrac- 
tor line. In generic models of potential-driven inflation, 
quantum fluctuations give rise to the phenomenon of eter- 
nal self-reproduction j(| 0, H, EJ E|. As soon as suit- 
able initial conditions occur in one domain, ultimately 
an infinite volume of space will undergo a long period of 
inflation followed by reheating. In this way, the eternal 
character of the self-reproduction process significantly al- 
leviates the problem of initial conditions for inflation. 

The main focus of this paper is to investigate self- 
reproduction in generic models of fc-inflation with La- 
grangians of the form We find that quantum "jumps" 
generically dominate over the "drift," and that self- 
reproduction is generically present. Self-reproduction in 
inflation can be studied using the stochastic ("diffusion") 
approach based on a Fokker-Planck (FP) equation (see 
Ref. 01 f° r a review). Although the FP approach de- 
pends on the choice of the spacelike foliation of space- 
time ("gauge"), one can use the FP equations to com- 
pute a number of gauge-independent characteristics, such 
as the fractal dimension of the inflating domain [llL 
and the probability of reaching the end-of-inflation point 
4> = 4>e of the configuration space. In particular, the 
presence of eternal self-reproduction is a gauge-invariant 
statement 0, . We show that the fractal dimension 
of the inflating domain is generically close to 3, confirm- 
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ing a strongly self-reproducing behavior in a wide class 
of models of fc-inflation. 

We also investigate the "stationarity" of the distribu- 
tion of 0, in the sense of Refs. There exists 
a (model-dependent) boundary value max such that the 
stochastic formalism breaks down for > 0max- Impos- 
ing a boundary condition at <f> = max , we find that the 
stationary volume-weighted distribution of the field i s 
concentrated near max . In the terminology of Ref. 
this constitutes a "runaway diffusion" behavior, mean- 
ing that the volume-averaged value of the field is driven 
to the imposed upper limit. Finally, we follow the evo- 
lution of the field along a single comoving worldlinc, 
starting at a generic value 0o < max - We find that the 
field will either reach the boundary — max or, with 
roughly equal probability, exit the inflationary regime 
through the end-of-inflation (reheating) point = 4>e- 
An exit into reheating becomes overwhelmingly proba- 
ble only if 0o < 4>r, where 4>r <C max is a model- 
dependent threshold value. Therefore, the semiclassical 
description of the inflationary evolution along a single co- 
moving worldline will remain valid only if the initial value 
0o is chosen within the interval 0# < 0o < 4>R- It is in- 
teresting to note that a similar conclusion was reached 
in a study of stochastic back-reaction effects in models 
of quintessence [Tfij . where the initial value of the infla- 
ton and the quintessence was found to be constrained to 
certain intervals. We compare this situation with the be- 
havior of models of potential-driven inflation, where such 
a threshold 4>r <SC max is absent and the exit into reheat- 
ing has probability 1 for all 0o ^ 0m ax , where max is the 
Planck boundary. 

Thus, the present investigation confirms that generic 
models of fc-inflation exhibit standard features of eter- 
nal inflation, and shows that strong effects of self- 
reproduction constrain the choice of initial conditions. 

II. MAIN RESULTS 

In this section we present the main line of arguments 
and the principal results of this paper. Full details of the 
relevant derivations are delegated to sections IlIII to Ivl 

In the cosmological context, we consider a spatially 
homogeneous field = 4>{t) and a flat FRW metric, 

g^dx^dx" = dt 2 - a 2 (t)dx 2 , (3) 

where a(t) is the scale factor. We also assume that the 
energy density of the field 0, 

e(X,<p) = 2X PiX -p, (4) 

dominates all the other forms of energy in the universe. 
(Partial derivatives are denoted by a comma, so p } x = 
dp/dX.) Then, as is well known, the Einstein equations 
yield 

d r 2 8?r ,r\ 

— — tl = KWe, K = — — = - ; „ , (5) 
a ' 3 3M 2 ,' v ' 



while the equation of motion for the field 0, 

— (p jV ) + ZHp^ v = p,,/,, v = 4> = V2X, (6) 

can be rewritten as 

e = -3H(e + p). (7) 

In Sections IIII AlIIlTJl we determine the conditions for 
the existence of attractor solutions in general models of 
fc-inflation. In a typical scenario, the field starts at a 
large initial value = 0o with < and gradually ap- 
proaches a smaller value = 4>e where inflation ends. 
If an attractor exists, the trajectory in the phase space 
quickly approaches the attractor curve; the equation of 
motion for 0(f) is effectively reduced to a first-order equa- 
tion. A slow-roll parameter can be introduced to quantify 
the deviation of the spacetime from exact de Sitter. In 
Sec. IIII El we show that a cosmological model with a La- 
grangian of the form admits an inflationary attractor 
solution under the following (sufficient) conditions: 

(i) The function K(<f>) remains positive at — > oo and 
is such that f°° -J K((f>)d(j) = oo (the integral diverges) 
and lim^oo^ ~ 3 / 2 K') = 0. For instance, functions 
K(<p) ~ s with s > —2 satisfy these two conditions. 2 

(ii) The equation L'(X) — has a unique root X > 
such that the energy density, s(Xq, 0) = —K(<j))L(Xo), is 
positive and bounded away from zero for large 0. (This 
condition excludes K((f>) ~ s with s < 0.) The attractor 
solution is X — X*(0), where the function A*(0) tends 
to Xo for large 0. 

(iii) The speed of sound of field perturbations is real- 
valued, c 2 > 0, for X — X and large (i.e. perturbations 
on the attractor background do not cause instabilities). 
It will be shown that c 2 > if L"{X )K'(</>) > for large 
0. The speed of sound c s (which is small in typical models 
of fc-inflation) plays the role of a slow-roll parameter. 

The attractor solution will be = v*(0); as long 
as the system is on the attractor, all the cosmological 
functions can be expressed as functions of 0. We ob- 
tain a complete asymptotic expansion of u*(0) at large 
(Sec. IIII Rll and a representation in terms of an inte- 
gral equation (Sec. IIII C|) . We show that the slow-roll 
condition generally holds (Sec. IIII F|) . Explicit asymp- 
totics in the leading slow-roll approximation are given in 

secnnm 

In the stochastic approach to self-reproduction, one is 
interested in averaged values of the field on sufficiently 
large distance scales, on which the fluctuations freeze and 
become uncorrelated. We show in Sec. II V Bl that the 



2 Technically, the two conditions arc independent, as exemplified 
by the functions K(<j>) = tj)' 1 and K{<j>) = (</>ln <j>)~ 2 . Note that 
the borderline case K(<f>) = 4>~ 2 must be considered separately 
(the attractor has the exact form tf> = const). Below we shall 
only consider K(<j>) ~ </> 3 with s > 0. 
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relevant scale L is not the Hubble horizon H 1 but the 
sound horizon, L ~ c s P _1 , where the speed of sound is 



V £ ,x 



(8) 



Let us denote the averaging on these scales by an overbar. 
The dynamics of the averaged field is a superposition of 
the deterministic motion and quantum "jumps" due to the 
backreaction of quantum fluctuations exiting the (sound) 
horizon. Since the attractor solution is approached expo- 
nentially quickly by other solutions, any deviations from 
the attractor are quickly suppressed. So we may assume 
that at every place the cosmological quantities such as H, 
c s , are always given by the functions of (b derived from 
the attractor solution fScc. IIII 0*|l . 

The backreaction of quantum fluctuations on the av- 
eraged field <j> may lead to a qualitative change in the 
evolution of the field, compared with the deterministic 
trajectory <j> = V*(<f>). To describe the dynamics of the 
averaged field during inflation, we use the stochastic or 
"diffusion" approach (see Ref. [T(| for a review) . The dis- 
tribution P c {4>, t) of values of <j> per coordinate (or comov- 
ing) volume satisfies the Fokker-Planck (FP) equation, 



dPc 
dt 



d 2 



(DPc) 



d 



(«,P C ) = L C P, 



(9) 



where u*(0) is the deterministic attractor trajectory and 
D{4>) is the "diffusion coefficient" describing the magni- 
tude of fluctuations. The computation of D{<f>) is per- 
formed in Sec. IIV CI and the result (in the leading slow- 
roll approximation) is expressed as 



D{4>) w const 



(10) 



The comoving distribution P c can be interpreted as 
the probability density, P c (<p,t)d<f>, for the value of <p at 
time t along a randomly chosen comoving worldline x = 
const. Also of interest is the distribution of <f> weighted by 
proper 3-volume within 3-surfaces of equal time t. This 
volume- weighted distribution, denoted P p ((f>, t), satisfies 
an analogous FP equation, which differs from Eq. (jj|J only 
by the presence of the expansion term 3P(0)P, 



BP d 2 

- — (DP 

dt ~ deb 2 1 p 



d 

— (u.Pp) + 3HP P = L p P. (11) 



These FP equations are supplemented by boundary con- 
ditions. One imposes an "exit-only" condition at the 
end-of-inflation point, <p = 4>e, and either a reflecting 
or an absorbing boundary condition at a suitable value 
4> = ^max (see Sec. IV Cl for explicit details). One then 
expects that solutions of Eq. J3J should decay with time, 
indicating that at late times a given comoving worldline 
will have almost surely exited the inflationary regime. In 
Sec. IV CI we show that the differential operator L c in 
the r.h.s. of Eq. 10 indeed has only negative eigenval- 
ues, indicating a late-time behavior P c (</>,<) cx e _At with 



A > 0. On the other hand, solutions of Eq. (|llfl may 
either decay or grow at late times, depending on whether 
the largest eigenvalue A max of the operator L p is negative 
or positive. A positive value of A max indicates the pres- 
ence of eternal self- reproduction (eternal inflation), with 
a stationary distribution P p ((b,t) cx e Amax *. It is found in 
Sec. IV Bl that A max > in generic models of /c-inflation, 
and that the distribution P p (<fi,t) is concentrated near 
the upper boundary <fi = max . Moreover, we show that 
A max ~ 3P max , where £f max = H(<j> max ) is the largest 
value of H(<j>) within the allowed range of <b. This shows 
that the inflating domain grows almost at the rate of the 
Hubble expansion. It is known that the inflating domain 
can be visualized as a self-similar fractal set [TTlll2fl . The 



fact that An 



3P max means that the fractal dimension 



of the inflating domain is close to 3, indicating a regime 
of strong self-reproduction. 

Finally, we calculate the probability of reaching the 
end-of-inflation boundary (j>E if the evolution starts at 
an initial value <f> — <j)Q (Sec. IVDjl . In all models of k- 
inflation, the consistency of the diffusion approach re- 
quires one to limit the allowed range of <fi by an upper 
boundary ^ max , even if the energy density does not reach 
the Planck scale fSec. HVE|) . The boundary cf) max is anal- 
ogous to the Planck boundary in models of potential- 
driven chaotic inflation, in that the semiclassical ("dif- 
fusion") approach fails for <f> > max because fluctua- 
tions become too large. One expects that the probability 
of exiting through the end-of-inflation point <Pe, rather 
than through the upper boundary <^ max , should be almost 
equal to 1 , indicating that almost all comoving observers 
will exit inflation normally. We show that there exists a 
threshold value 4>r <C </> max such that the exit through 
4>e is overwhelmingly likely for initial values 0o < 4>r- 
(Such a threshold is absent in potential-driven inflation.) 
Hence, a given comoving worldline with 4>q < 4>r will al- 
most surely (with probability 1) eventually exit inflation 
and enter the reheating phase. Nevertheless, the total 
proper 3-volume of the inflating domains will grow expo- 
nentially with time. Thus, we demonstrate that generic 
models of /c-inflation exhibit the standard features of eter- 
nal self-reproduction. 



III. ATTRACTOR SOLUTIONS IN 
fc-INFLATION 

In this section we study the cosmological dynamics of 
the spatially homogeneous field <p(t) with a Lagrangian 
p(<j>, (b) and investigate the existence of attractor solu- 
tions in the regime <j) — > oo. We show that attractor 
solutions exist in a broad class of models and determine 
the asymptotic forms of the attractors. As a particular 
example, we select the Lagrangian 



p(0,0) = X(0)Q(0). 



(12) 



Since the equation of motion © does not depend ex- 
plicitly on time t, it is convenient to consider solutions 
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(j)(t) as curves in the phase plane (0, v), where we denote 
v = <p. The equation of motion in the phase plane has 
the form 



-j7 = g{v,<t>), 



(13) 



where g(v,4>) is an auxiliary function expressed through 
the Lagrangian p(v, <p) as 

g{v,(f) = — f 3k yj e(v, <j>)p iV + vp. v<t> - pA , (14) 

vp, vv v ' 

e(v,<j>) = vp jV -p. (15) 

In the particular case of the Lagrangian H12|) , the function 
g(v, (/)) is given by 



vQ"(v) vQ"(v) K{(p) 



(16) 
(17) 



Here and below we assume that K((f>) > when 4> — » oo. 
We can always change the sign of Q(0) if this is not the 
case. Note that K(<f>) may have a root at a finite <j>, 
e.g. at the endpoint of inflation, but we only consider the 
inflationary regime where <p is sufficiently large and K(4>) 
remains positive. 
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Figure 1: A numerically obtained phase plot for the La- 
grangian p(X,(j>) — (—X + X )arctan0. Trajectories start- 
ing at large negative <j> with large positive v = (j> approach an 
attractor solution. The attractor trajectory approximately 
coincides with the line v = vo for sufficiently large \<f>\. 



p(v, 4>) is exponential in v. Hence, the growth of a non- 
attractor solution for large <f> is at least exponential, both 
for small and for large deviations from the attractor. 

For convenience, below we always consider the quad- 
rant 4> > and v < in the phase plane, since the 
opposite choice cf> — > — oo and v > is treated completely 
analogously. In the general case, we expect that the at- 
tractor solution u*(</>) will have a relatively slow behav- 
ior at 4> — > oo, compared with an exponential behavior 
of nearby non-attractor solutions. Thus, a suitable con- 
dition for selecting the attractor solution is that u*(0) 
should grow slower than exponentially as — > oo, while 
nearby solutions grow exponentially. A formal way to 
express this condition is 



Definition of attractors 



lim 



lnv*(tf>) = 0. 



(18) 



Let us now motivate the definition of an attractor by 
analyzing the typical behavior of solutions in terms of 
the function v(<f). A sample set of trajectories in the 
phase plane (4>, v) is shown in Fig.^ Trajectories starting 
at large (f> and v < will quickly approach an almost 
horizontal line, v ~ vq — const, and then proceed along 
that line more slowly towards <p = 0. It is intuitively 
clear from the figure that the line v ~ vq (rather than 
any of the neighbor trajectories) should be considered 
the attractor solution. To make this statement precise, 
we need a formal criterion that would distinguish the 
attractor line from nearby solutions. 

We begin by presenting a heuristic motivation for this 
criterion (which will be Eq. I|18|) below). In Fig. ^ 
the attractor solution stays approximately constant as 
4> — * — oo, while a generic neighbor trajectory vq + Sv(<p) 
quickly moves away from the attractor. Let us there- 
fore examine the growth of 6v(<f)) with 4> as <f> — > — oo. 
If v w vq is a solution of Eq. I|13fl . then g(vo,(f>) ~ 
and g(vo + Sv, 4>) ~ 9,v(vo, 4>)6v. Since -j^Sv w g. v Sv, an 
initially small deviation Sv((f>) grows exponentially with 
4> if g,v(vo,(j)) 0. Assuming a power-law dependence 
p(v, 4>) ~ v n for large v, one finds from Eq. I|14|) that 
the function g(v, (j)) will grow at least linearly with v at 
fixed <f>; the growth of g(v, <fi) with v will be even faster if 



We shall use Eqs. I|13|) and (|18|l as the definition of the 
attractor solution w*(0), together with the assumption 
that nearby trajectories do not satisfy Eq. (|18fl. We stress 
that the attractor solution is singled out by an asymptotic 
condition at <p — > oo. At a finite value of cf>, all the 
trajectories approach each other and no single solution 
appears to be special. 

Let us remark that the definition (|18|l may be unduly 
restrictive. In particular, it does not allow attractors that 
are approached more slowly than exponentially. In such 
cases, the definition H18|) may still be used to select the 
attractor solution after a suitable change of variable <fi — > 
(j). Moreover, it may happen that either every solution or 
no solution of Eq. (|13fl satisfies Eq. (|18fl . In such cases, 
a more detailed analysis is required to investigate the 
existence of attractors. Nevertheless, the condition i|18|) 
is adequate for the analysis of /c-inflation. 



B. Asymptotics in the general case 

The next task is to determine the possible asymptotic 
behavior of the attractor «*(</>), given the equation of 
motion i|13|) . Let us assume that the function g[y,4>) has 



5 



an asymptotic expansion at 



oo of the form 



g(v,<f>)=Ao(v)Bo{<l>)+Mv) B i(4>) 



(19) 



where -Bi(</>) is subdominant to Bq(<P) as (p ~ > oo. The 
asymptotic expansion of g{y, (ft) at cj> — > oo will have the 
form l|19l) if the Lagrangian p(y, 4>) is polynomial in v. In 
particular, Eq. Ijltil) is manifestly in this form. 

In most cases, the term Ao(v)Bo((j)) in Eq. (fTT)|l will 
dominate at large <p, so that an approximate general 
(i.e. non-attractor) solution v{4>) of Eq. i|T5|) is found as 



dv' 



(20) 



where v(cj>i) = v% is an arbitrary initial condition. The 
technical assumptions for the dominance of the term 
Ao(v)Bo(cf>) are that both Ao(v) and Ai{v) are bounded 
for finite v and that Ai(v) does not grow faster than 
Aq(v) for large v. It is straightforward to verify, using 
Eq. Q14JI. that these assumptions hold for a Lagrangian 
p(<f>,v) which is polynomial in v. 

To study the behavior of the general solution (|2U|l for 
large </>, we need to consider two mutually exclusive cases: 
either J B (4>)d(j) converges at <p —* °°, or it diverges. 
If this integral converges, there exist solutions v(4>) that 
approach any given value v\ as cj> — > oo, namely, 



Vl dv' 
(<t>) Aq(v') 



(21) 



Since v% can be chosen continuously, there are infinitely 
many solutions v(t) that tend to a constant limit at 
<t> — > oo. It follows that there is no unique attractor solu- 
tion in this case. In the second case, f°° Bo((f>)d(f) = oo, 
there may exist solutions v{cj>) approaching a root of 
Aq(v) as 4> — ► oo. To simplify the analysis, let us change 
the variables, 



Bo 



(22) 



noting that tf> — ► oo together with — * oo. We can now 
rewrite the equation of motion through the new variable 
<p and look for attractor solutions at (f> — > oo using the 
same criteria as before. The asymptotic form of g(v, <p) = 
g(v,cj ) )/B (cb) will be 



){v,4>) = A {v) + Ai{v) 



B (4>) 



(23) 



Essentially, the replacement <f> — > <p is equivalent to the 
assumption that Bq{4>) = 1 in Eq. (|19fl . This is what we 
shall assume in this and the following sections. 

If v*((j)) is an attractor, Eq. 1)18(1 with Bo(<p) = 1 gives 



= lim 



v*(<f>) 



lim 

d> — >oo 



(24) 



and hence v*{4>) approaches a root of A (v)/v as <f> — > 
oo. Thus, roots of the function Aq(v)/v correspond 
to solutions that are approximately constant for large 
(j>. Our previous analysis shows that such solutions, 
v{4>) ~ v , will be attractors if A' (vo) > and repul- 
sors if A' (vo) < 0. Therefore, we focus attention on the 
case when the function Aq(v) has a finite root vq < 
such that A' q (vq) > 0, and we look for an attractor so- 
lution that approaches v = vq when — > oo. 3 
Additionally, a physically meaningful inflationary attrac- 
tor solution v*(0) must have a positive energy density, 
e (v *((/>), <f>) > 0. Thus the root must also satisfy the 
condition 



lim e(vo, <ft) > 0. 



(25) 



Since V*((f>) « const for large <p, an approximation to 
the attractor solution can be found heuristically by ne- 
glecting v' in Eq. (|13JI . Clearly, this is equivalent to the 
frequently used slow-roll approximation where one ne- 
glects <p and keeps only <p in the equations of motion for 
4>(t). Thus, the slow-roll approximation v sr {<f) to the at- 
tractor solution is the root of g(v sr (4>), <j)) = which is 
near vq as ^ oo. 

We can now obtain a complete asymptotic expan- 
sion of the attractor solution «*(</>) at <f> — > oo. Since 
lim^oo g tV (y, 4>) = A' (v) ^ in a neighborhood of vq, 
we may solve Eq. ()13|l with respect to v at least in some 
neighborhood of (cj) = oo, v = vq) and obtain 



(26) 



where g -1 (-,0) is the inverse of g(v,4>) with respect 
to v. (If the equation g{v,<p) = has several roots 
with respect to v, an appropriate choice of the branch 
of the inverse function must be made to ensure that 
g~ (0, 4>) — * vq as 4> —>■ oo). Let us apply the method 
of iteration to Eq. I|2t)|) , starting with the initial approx- 
imation f(o)(0) = const. The first approximation will 
coincide with the slow-roll solution, 



v (i)(4>) = 9 X ($>A) = v sr {(f)- 



(27) 



3 Thus, we omit the consideration of cases where the existence 
of attractors is less straightforward to establish. For instance, 
the function Ao(v)/v may have no finite roots, may approach 
zero as » — > oo, may have a root vq = 0, or a higher-order 
(multiple) root. Also, an attractor solution may exist not due to 
a root of Ao(v) but due to a cancellation among several terms 
in the expansion 1191 . Finally, an attractor solution may grow 
without bound, v *(</>) — » oo as c/> — * oo, instead of approaching a 
constant. In every such case, a more detailed analysis is needed 
to investigate the existence of attractors, possibly involving a 
further change of variable, such as v((f>) = v(<j>)F(<f>), where F(<p) 
is a function with a suitable asymptotic behavior at <j> — > oo. See, 
for instance, AppcndixlAlwhere we use this technique to establish 
the existence of unique attractors in models of potential-driven 
inflation. 
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We shall now show that the successive approximations 
v (n)(4')i n — 1>2,..., differ from the exact attractor by 
terms of progressively higher asymptotic order at <fi — > oo. 
It will follow that the sequence {v n (<p)} converges to the 
attractor in the asymptotic sense (i.e. for <j) — ► oo at fixed 
n). This asymptotic approximation is similar in spirit to 
the asymptotic expansion using a hierarchy of slow-roll 
parameters, which was developed in Ref. 17). 

By construction, successive approximations vr n \(<fi) 
satisfy 



dV{ n ){4>) 



(28) 



We denote by (<j>) the exact attractor solution satisfying 
Eq. l(2l)|) . Then the deviation of U(„) from v» satisfies the 
recurrence relation 



V(n+1) ~ V* 
1 



gA v *A) 



■g-^v'^^-g- 1 ^) 

[v(n)-V*],+0 (v {n) - 



(29) 



Since v^(<p) — v*(4>) decays at (f> — > oo (this is so be- 
cause both functions approach the same constant vo), it 
follows that (i>(„) — v») also decays at <f> — > oo for all 
n. Hence, the terms («(„) — w*),0 and («(„) — u») 2 have a 
higher asymptotic order (i.e. decay faster as (f> — > oo) than 
( v (n) ~ v *)- Since for large <f> we have g, v {v*) ~ A' (vq), 
which is a nonzero constant, it follows that (w( n+1 ) — v*) 
has a higher asymptotic order than (i>(„) — U*) at 4> — > oo. 



C. Existence of attractors 

In the previous section, we have shown that an asymp- 
totic expansion of the attractor solution may be obtained 
by iterating Eq. H2(jf> . However, the successive terms of 
the expansion will contain derivatives of Bi(<p) of suc- 
cessively higher orders. Since the n-th derivative of an 
analytic function usually grows as quickly as ~ n!, the 
asymptotic expansion could be a divergent series. Such 
a series cannot be used to reconstruct the attractor so- 
lution exactly and establish its existence. In fact, the 
asymptotic property of the sequence i!( n ) (</>) was derived 
assuming that the exact attractor solution u*(</>) exists. 

The existence of the attractor solution can be estab- 
lished by producing a convergent sequence based on an 
integral equation. To this end, we rewrite the function 
g(v,<j>) as 

g(v, <t>) = {v- v ) a+(v- v Q ) 2 A(v) + B(v, 0), (30) 

where a = A' (vo) > and A(v) , B (v , (j>) are auxiliary 
functions such that A(y) is regular at v = vq and B(v, 4>) 
decays as <p ~~ * 00 uniformly for all v near vq. (It is 
possible to express g(v,4>) m the above form because 
g,v(vo,4>) — > ol as cj) — > oo. The function B(v,4>) has 
the stated property since A\(v) is regular at vq.) A triv- 
ial redefinition — > a~ 1 (f> effectively sets a = 1, which 



simplifies the analysis. Then Eq. (|13[1 with the bound- 
ary condition v(oo) = vq is equivalent to the following 
integral equation, 

poo 

v{4>) =v Q -e* e~* [ {v{4>') - v f A{v{<p')) 
J 4, 

+ B(v(<p'), <t>')] (31) 

This equation can be iterated starting from the initial 
approximation v^{(j>) = vq, which produces successive 
approximations t)W(^), n = 1,2,... We shall now show 
that the sequence v^ n ' (4>) converges as n — > oo for suffi- 
ciently large (but fixed) cj). The limit of the sequence will 
be, evidently, the attractor solution i>* (</>). 

We shall prove the convergence of the sequence 
{ v( n \cj))} by showing that the difference between succes- 
sive approximations, |u( n + x ) ((f)) — u' n )(0)|, decays with n 
at fixed <j). Let us first estimate the difference A n (<fi) = 
V^ n \cj)) — vq at some fixed cj). By construction, we have 

A„+i(0) 

p CO 

= -e* / e~*' A^( w (")(0')) + S(t-( n )(0'),0') #'• 

(32) 

To determine a bound for A„ using this equation, we need 
some bounds on the functions A(v) and B(v, <f>) valid near 
v — vq and for large <fi. By assumption, B(v,(j)) decays 
uniformly in <p, so for any small number 8 > we can 
choose large enough, say cj) > </>o, so that 



\B(v,<j))\<8 
\B iV (v,<f>)\ < 8 



for \v — vo\ < 8 and <f> > (j) . (33) 



Since A(v) is regular at v — vq, there exists a constant 
C > such that \A(v)\ < C for \v - v \ < 8. Then the 
difference A n+ i is estimated using Eq. (|32|) as follows, 



|A„+i(^)| <C|An(0)r + d 

Since the initial approximation is v^ ((j>) 
to see that, with 8 < (4C) _1 , we have 



(34) 

vq, it is easy 



|A»(0)| < 



1 - VI - 4C8 
2C 



< 28, for Vn and <fi > 4>q. 



(35) 

Let us choose 8 < (4C) 1 and find a corresponding value 
of 4>q such that Eq. (|33|l holds. Then we compute 



,(«+!) 



(0) 



,(n) 



(0) 



< 



c 



v (n)_ v (n-l) 



(36) 



Since 



- t-o) 2 - [v {n - 1] - V Q f 



(„(») _ - Vo + - Uo ) 



< 4 „(«) _ S. 



(37) 
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we have 



< 5(5 



(38) 



Now we use induction in n, starting with 



< 



< 



e*-* '6 = 5, 



(39) 



to establish the bound 



- < (55)" (5, for Vn and > <fo. 

. (40) 

This bound tends to zero for n — > oo as long as 5 < g. 
If necessary, we can diminish <5 so that 5 < g and choose 
a corresponding value of (f>o such that Eq. (|33() holds. 
Then the sequence v^ n \4>) converges as n — > oo at any 
fixed 4> > 4>q. 



D. Summary of assumptions 

In the preceding sections we have established that a 
general model of fc-inflation admits an attractor solution 
v*(4>) under the following (sufficient) conditions: 

1) The function g(v, <j>) defined by Eq. (|14f> has an 
asymptotic expansion at <fi — > oo of the form (|19(l , where 
Bi((j}) I Bq((/)) — ► as cj> — > oo . This asymptotic expansion 
defines the functions A (v), Ai(v). 

2) The function Aq(v) has a unique root vq < such 
that A' (v a ) > 0. 

3) The function A\(v) is regular at vq and does not 
grow faster than Aq(v) for large v. 

4) The function Bo(<p) is such that J Bo((f>)d(f) = oo 
(the integral diverges). 

We have shown that the attractor solution v»(0) ap- 
proaches vq as <j) — > oo and has a well-defined asymptotic 
expansion for (f> — > oo, which can be determined by iter- 
ating Eq. H26fl starting with w(0) = Uq- The first iteration 
yields the slow-roll approximation to the attractor, which 
is the function v sr ((f>) defined by 



g{v sr {<j>), <j>) = 0. 



(41) 



Assuming that (generically) Ai(vq) ^ 0, we can obtain 
the following leading-order approximate expression for 
the slow-roll solution, 



V sr (4>) « Vq 



A' Q {v ) B {<t>Y 



(42) 



In addition to the above assumptions, the following 
two conditions must be satisfied in order for the attractor 
solution to be physically relevant: (i) The energy density, 



e(v,(p) = vp >v — p, is positive for all v near v Q and does 
not approach zero as 4> — > oo. (ii) The speed of sound is 
real-valued, c 2 s > 0. The condition (i) will be satisfied if 



e(vo, <f>) > for all sufficiently large • 
says that 



The condition (ii) 



P.v 



P,v 



VP,; 



>0. 



(43) 



An explicit form of the condition l|43|) for the La- 
grangian l|12|l is given in Eq. (|52|l below. 



E. Existence of attractors in fc-inflation 

We now apply the preceding constructions to a model 
with the Lagrangian i|12|) . We first need to analyze the 
asymptotic form of the equation of motion, v' = g(v, <j>), 
at <j) — > oo. The function g(v,<j>) is given by Eq. (|16f) . 
which is in the form (|19fl with only two terms, 

g{v,4>) = A Q {v)B {4>) + A 1 (v)B 1 ( ( f>), (44) 

A (v) = -3k v , A 1 (v) = — (45) 

vQ"(v) vQ"(v) 



i(v) = vQ'(v)-Q(v); 



(46) 



B {4>)= y/K{$), B 1 (0) = — [In K{4>)] . (47) 

The asymptotic behavior of the functions Bq(4>) and 
B\ ((f)) at 4> —> oo may belong to one of the three cases: (i) 
B (<t>) « Bx{4,), (iij B (<t>) » Bx{f), and (hi) B {4>) - 
Bi((f>). In the first case, the equation of motion is dom- 
inated by A\(v) which cannot have roots Ai(vq) = 0, 
otherwise the energy density would become small on the 
attractor, e(t;*,0) w s(vq,4>) = at 4> ~~ * °°j an d could 
not dominate the energy density of other matter. Hence, 
in case (i) inflationary attractors do not exist. Case (ii) 
covers all functions K(cf>) that either grow with (f>, or tend 
to a constant, or decay slower than <fi~ 2 as <f> — > oo. Case 
(iii) corresponds to a specific family of Lagrangians with 
K(4>) oc <f>~ 2 , which gives B ((j)) ~ -Bi(</>) oc </> _1 . After 
a change of variable, cf> — > <f> = In cf), this case is reduced 
to case (ii) with Bq = 1 and B\ = 0. Therefore, we shall 
confine our attention to case (ii). The precise conditions 
for that case are 



Bi(4>) _ K'(4>) 



B (cb) K^icji) 



as 



Bo 



(48) 
(49) 



Assuming a polynomial growth of Q(v) ~ v n for large 
v, where n > 2, we find that A$(v) grows faster than 
Ai(v); this verifies condition 3 in Sec. IIII Dl Therefore, 
the term Aq(v)Bq((/)) dominates the equation of motion 
for non-attractor solutions at large 4> and, as shown in 
previous sections, an attractor solution v*(4>) must tend 
to a root vq 7^ of Aq(v) as 4> — > oo. The function Aq(v) 
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is expressed by Eq. I|45|l . Since, as we already found, 
e(v) must be nonzero to drive inflation, it follows that v 
must be a root of Q'(v). Let us assume that v < is 
such a root, Q'(vq) = 0. It remains to verify the other 
conditions for the existence of the attractor. 

The energy density is positive on the attractor solution, 
c = -K((t>)Q(v ) > 0, if K{$) > and Q(v ) < 0. 
Thus, the only admissible roots of Q'(vo) are such that 
K(4>)Q(vq) < 0. The condition ^(^o) > is satisfied 
since v < 0, e > 0, and thus 



Q"(vo) 
v Q"(v ) 



-3k- 



> 0, (50) 



where we have assumed that (generically) Q"{vq) ^ 
0. Finally, the speed of sound must satisfy the condi- 
tion (O, 



P.r 



"0 



> 0. 



(51) 



This condition effectively constrains the deviation of the 
attractor solution v*((j>) from v — v = const, requiring 
that v*(t) < vq- Because of that, we have Q'(u*) < and 
g(v*,(f>) > 0. Then the equation of motion, dv*/d(j) = 
g(v*,4>), entails that the negative value t>*(</>) monoton- 
ically increases (in the algebraic sense) as (j> decreases, 
at least until the slow-roll approximation v*(cj>) w vq 
breaks down. Numerical simulations for some model La- 
grangians show that v* (4>) is typically monotonic in 4> all 
the way until the end of inflation. 

Using the approximate solution (|42[1 and the fact that 
vq < 0, we rewrite the condition i|51|) as 



1 A^vo) B^) = Vf 

\v \ A' (v ) B (<f>) 3k \v \ Q"(v ) K»/*{4>) 



> 0. 
(52) 

It follows that c 2 s > if Q"(v )K'(<j>) > for large <j>. 
We conclude that a model of fc-inflation satisfying the 
conditions listed in Sec.[n]will admit an attractor solution 
with properties suitable for inflationary cosmology. 



F. Slow-roll condition 

In this section we show that attractor solutions in k- 
inflation models satisfy the slow-roll condition for large 
<p. The slow-roll condition means that the change in the 
Hubble rate, AH = HAt, is negligible during one Hubble 
time At = H~ l , i.e. 



H 



< H. 



(53) 



This condition guarantees that the spacetime is locally 
sufficiently close to de Sitter during inflation. 
It follows from Eq. Q that 



1 dH 
H~dt 



1 de 
e dt 



3 -H 
2 



e + p 



(54) 



therefore the slow-roll condition with the Lagrangian i|12|) 
is equivalent to |e +p\ <C £ or 



Vp,v 




vQ'(v) 


VP.v - P 




vQ'(v) - Q{v) 



< 1. 



(55) 



Since the attractor approaches a root vo of Q'(v) as 
cj> — > oo, and since Q[vq) < by the condition of the 
positivity of the energy density, the slow-roll condition 
will be satisfied for large enough <f>. More precisely, <fi 
must be large enough so that 



v Q Q"(v ) 
Q(vo) 



\v(<f>)-va\ < 1. 



(56) 



The slow-roll condition is violated (and inflation ends) 
when d — 0, which is equivalent to H + H 2 = or e + 
3p = 0. This occurs at a time t when either vQ'(v) + 
2Q(v) = with v = <j>(t), or K(<f>(t)) = (then e = 
p = at the same time). Generically, we may expect 
that the former condition does not take place; so the 
end-of-inflation point is the root 4>e of K(4>e) = 0. For 
instance, this is the case when K((f>) is monotonic and 



Q(v) 



-C\V 



-c 2 v n2 with ci > 0, C2 > 0, n% > n\ > 1. 



G. Asymptotics of attractors in fc-inflation 

Assuming that the conditions for the existence of the 
attractor are met, let us now compute the asymptotic 
form of the attractor solution for large 4>. All the func- 
tions describing the homogeneous cosmological solutions 
(</>, H, c s , etc.) can be expressed as functions of <p. Sub- 
stituting Eqs. P5 jl -l(i7 jl into Eq. we find 



4> = v*(<f>) w v - «i (</>), 



3kQ"(v ) A-3/ 2 (0)' 
2 

0o + vot + 



-Q(vo) 



3kQ"(v ) A/ K(cb Q + v Q t)' 



H(<f>) « ^~K{ct,)Q{v,) (l + ^M vM 



(57) 
(58) 

(59) 
(60) 



In a{4>) 



H 



v 



<t>0 



«0 



VKtfjty, (61) 

v^QW K'{<j>) 



v Zk\v \Q"{v )K3/ 2 ^Y 

The assumed conditions are 

v < 0, Q> ) = 0, Q(v ) < 0, 
K(4>) > 0, K'(<f>)Q"(v ) > 0. 



Note that c s «l for large 4>, due to the condition 
The small value of c 2 a can be considered a "slow-roll pa- 
rameter," i.e. a parameter describing the smallness of the 



(62) 



(63) 
(64) 
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deviation from the exact de Sitter evolution. The slow- 
roll condition (|56[l holds if 



VqQ"(vq) 
Q(vo) 



K>{cj>) 



< 1. 



(65) 



This inequality determines a model-dependent slow-roll 
range <f> > 4> sr . The slow-roll approximation becomes 
increasingly precise in the large <p limit where c s — > 0. 



IV. MAGNITUDE OF QUANTUM 
FLUCTUATIONS 



In this section we compute the magnitude of quantum 
fluctuations of the field <f>. This calculation will yield the 
relevant kinetic coefficients for the Fokker-Planck equa- 
tion, which will be considered in Sec. 



A. Quantization of fluctuations 

We consider cosmological perturbations of the field cf> 
on the inflationary attractor solution in a spatially flat 
FRW universe. We follow Ref. where the quantum 
theory of perturbations for fields with noncanonical ki- 
netic terms was developed. 

Perturbations of the field </> give rise to scalar pertur- 
bations of the metric. Their simultaneous dynamics is 
described by a single scalar field u with the (classical) 
equation of motion 



u" - cjAu u = 0, 

z 



(66) 



where c s is the speed of sound, the prime denotes deriva- 
tives with respect to the conformal time, 



du 



a(t) 



dr\ dt 
and the auxiliary function z(t) is defined by 



a V e + v 

c s H 



(67) 



(68) 



The fluctuation 8(f) of the field cf> and the Newtonian grav- 
itational potential $ are expressed through u (in the lon- 
gitudinal gauge) as 



4> 
H 



u <& 
~zH~ H' 

3/^ o . i d 

z 2 A~ 1 — 

2 dt 



(;)■ 



(69) 
(70) 



It is convenient to consider the Fourier modes 84>k, 
Uk of the perturbation variables; then the inverse Laplace 
operator in the above equation becomes simply — fc~ 2 . 



For a model with the Lagrangian l|12|) and the attractor 
solution found in Sec. 1111 Gl we have 



z 4'*Wii (1+0(a) ' ^ 



zH = \v \ a(t)^K(4,)QZ (1 + 0(c 2 s )) 



(72) 



where we denoted for brevity Qo = Q{vo), Q'o = Q"(vq)- 
We find that z(t) ss a(t) -const; here and below we denote 
by the approximation obtained by neglecting terms 
of order c 2 <C 1. Within this approximation, we can 
express the field fluctuation 54>k as 




It is known that the action functional for the field 
it(x, rj) has the canonical form. Therefore, one quantizes 
the field w(x, rf) by postulating the mode expansion 

— m —[a ic ut(rj)e i ^ + H.c.], (74) 
(2tt) V ^ 

where dk is a canonical annihilation operator satisfying 
[dk,a£.,] = 5(k — k'), and "i?.c." denotes Hermitian con- 
jugate terms. The mode functions Uk{rj) are normalized 
by the condition Im(ufcit£) = 1 and satisfy 



u'l + c 2 s k 2 u k - —u k = 0, 



a 
a 



(75) 



where we replaced z" / z by a" /a, which is justified within 
the slow-roll range. The conformal time variable rj can 
be expressed as a function of <fi, 



o Ma(</>) 



d<t> 
o I «o| 

1 



exp 



«VI<3o| 



Ky/\Q \K(<I>) 



exp 



Ky/\Qo\ 
\V \ 



It follows that a(rj) « {Hrf) 1 and a" /a 



-2 V - 



(76) 



so 



the mode functions Uk(r)) can be chosen as the standard 
Bunch-Davies mode functions for the massless field (ex- 
cept for the extra factors c s ), 



Uk{r)) 



1 



Jc 3 kr/ 



1 



c s krj 



(77) 



Since the field modes Sfik are linearly related to u^, 
a similar mode expansion holds for the quantum field 
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S<p{x,t), except for different mode functions as given by 
Eq. J7j|: 



[a k 5MvV k ' x + H.c] , (78) 



d 3 k 1 

(2tt) 3/2 71 

Uk 3/-ta 2 f 2 / Qo_ ^ ( u k 



ay/K(<l>)QZ 2fc 2 Y |Q | dt V a 



IQol 



3 vl£ 



Q'o 



Qg 2 1 - iCs fc7? V |Qo! 



(79) 



We may now disregard the term of order (? s and obtain 



'\Qo\ 



7?e 



1 



IQo 



^0 

(80) 

We conclude that in the slow-roll regime the mode func- 
tion 5(j>k(T]) is proportional to the standard Bunch-Davies 
mode function evaluated at the wavenumber ck. 



B. Scales of averaging 

We are now interested in the spatially averaged pertur- 
bations 5<fi of the field <j>. It is well known that quan- 
tum fluctuations of a scalar field (in the Bunch-Davies 
vacuum state) in de Sitter spacetime become uncorre- 
cted on Hubble distance and time scales, L ~ St ~ H^ 1 , 
and are transformed into classical perturbations. In the 
present case, the relevant distance scale is the sound hori- 
zon scale, L ~ c s i? _1 , while the time scale remains un- 
changed, dt ~ Although this was already noted 
e.g. in Ref. 0], we shall now justify this statement more 
formally using arguments similar to those of Ref. [19j . 

The averaged fluctuation field 5(f) can be defined as 

S4>{x,t)= /dV&^x'^W^L-^x-x'l), (81) 



where L is the averaging scale and W(q) is a suitable 
window function that quickly decays for q > 1. Then the 
effective noise field is defined as 



(82) 



For a minimally coupled, massless scalar field </>, it was 
shown in Ref. [ijj that correlators of the noise field £ have 
the following large-distance asymptotics, which hold for 
a wide class of window functions W: 



<£(x,iK(x',t')) 



-2H\t-t'\ 



H 4 |x - : 



(83) 



the field 6<j> differs from the minimally coupled mass- 
less field considered in Ref. 0] only by the replacement 
k — > c s k and by an overall factor; this is easy to see by 
examining the mode function (|8(jp. Therefore, the calcu- 
lations of Ref. 19] apply also to the present case if we 



rescale the length as L 



L. The time scale remains 



unchanged, whereas the distance scale is found from the 
relation c^ 1 L = H -1 , which yields L = c s -ff _1 . 



C. Diffusion coefficient 

The effective dynamics of the averaged field 
scribed by a Langevin equation, 



is de- 



(84) 



where £ (x, t) is the effective "noise" field, which encapsu- 
lates the effects of subhorizon quantum fluctuations. In 
the stochastic formalism, the "noise" is treated as a clas- 
sical field, i.e. a Gaussian random function with known 
correlations ("colored noise") 0, Eol EH ] . It was shown 
in Sec. IIVBI that the field £ is essentially uncorrelated 
beyond Hubble time scales St ~ H and beyond distance 
scales L ~ c s _ff _1 . To simplify the analysis, one may 
treat the "noise" field £ approximately as white noise as 
long as one considers variations of <f> only on sufficiently 
large distance and time scales. Such variations of </> can 
then be described by the following Langevin equation, 
which is this time a difference equation, 

0(x, t + Si)- 0(x, t) = (4>)St + ^2D(4>)StS~f, (85) 

where D is a "diffusion coefficient" and £7 is a fiducial 
Gaussian random variable with unit variance. The coef- 
ficient D{4>) is defined through the standard deviation of 
the quantum "jump" in the field <p during a time interval 
St, 



2DSt= <^ 2 (x,i)) 



L~c s H- 




(86) 

where we explicitly indicate averaging the field 4zS(/)(x, t) 
on distance scales L ~ c s i/ _1 . 

As long as one considers variations of <f> on scales at 
least ~ H^ 1 . the difference equation l|85|l is equivalent to 
the stochastic differential equation 



dt 



(87) 



where 7 is a fiducial "white noise" field. The above equa- 
tion is, in turn, equivalent to a FP equation for the time- 
dependent probability distribution P c (^, t), 



This indicates a quick decay on time and distance scales 
of order the Hubble horizon. In the present calculation, 
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where we use the Ito factor ordering in the FP equa- 
tion [13 . 

We can now compute the diffusion coefficient £>(</>) 
in the presently considered model of fc-inflation, using 
Eq. H86|) . Field fluctuations on distance scales L ~ c s iJ _1 
in the Bunch-Davies vacuum state are due to field modes 
exiting the sound horizon. Therefore, the magnitude of 
these fluctuations can be estimated by evaluating the ab- 
solute value of the mode function (|80|l at the wavenum- 
ber k that corresponds to the sound horizon crossing, 
c s k = aH w \rj\~~ . Since the quantum field -^S4> has the 
standard mode expansion with mode functions jjj6<fik, we 
find using Eq. (JSDJ) that 



D ic s krj 



c s ki] 



\Qo\ 

Qo 



'^TT7-^-^==ic s krj 2 H. 



(89) 



(90) 



Therefore, the typical fluctuation in -^5(j) averaged on 
scales L ~ c s _ff _1 and accumulated during one Hubble 
time, St = if -1 , is 



(54> 2 5t 2 ) 



L~c s H- 







dt Hk 









c s fcj r/\ — 1 



« 2 IQol 



(91) 



4^ 2 Qo s ' 

Finally, the coefficient D((f>) is computed from Eq. I|86l) 

k 2 \Q \H 



D 4 i 1 *'" 



Stt 2 Q' ' cf 



(92) 



Since c s <C 1 for most of the allowed range of 0, the 
diffusion coefficient is typically large. As we shall see, 
this leads to the presence of self-reproduction. 

A more explicit expression for D(<p) can be found using 
Eqs. (EOJ) and @£J|, 

^^^wTOBT, (93) 



8ir 2 



since c s (4>) — > 0, whereas -ff(0) ~ wK^j remains con- 
stant or grows as <fi ~ > oo. Therefore, diffusion will dom- 
inate over the deterministic drift for all sufficiently large 
cj). Moreover, since the slow-roll condition is essentially 
the same as c 2 <C 1, diffusion becomes negligible only 
near the end of the slow-roll range or even beyond that, 
near the end of inflation. A wide range of </> where dif- 
fusion dominates over the deterministic drift is a quali- 
tative indication of eternal self-reproduction. A precise 
criterion for the presence of self-reproduction is provided 
by the Fokker-Planck equation for the volume-weighted 
distribution of the field, which will be analyzed in Sec.l\H 



E. Validity of diffusion approximation 

We have computed the magnitude 8<p of a typical ran- 
dom "jump", defined as 5<fi ~ y / 2D((f))5t, during one Hub- 
ble time St = H^ 1 . The diffusion approximation is valid 
if 5<j) is "small" in the sense that it causes small changes 
in the energy density of the field, 



5e = e^Scf) < e(4>). 



(96) 



If this condition does not hold, we would not be justified 
in assuming that the "jumps" cause a perturbatively small 
change in the local expansion rate H. Since s^/e ~ 
H^/H, an equivalent condition is 



IT 



^y/2D(J>)H-i « 1. 



(97) 



There is another aspect that limits the validity of the 
diffusion approach. As we have seen, the "jumps" are un- 
corrected on distances of order c s iJ _1 , which is signifi- 
cantly smaller than the Hubble scale H^ 1 since c s « 1. 
A comoving scale c s iJ _1 grows to Hubble size within 
the time 5t g = H^ 1 IncJ 1 . Thus, a single Hubble-sized 
region will contain small inhomogeneities in the values 
of <f> on scales above c s -ff _1 , generated during the time 
interval 8t g > H~ l . The typical magnitude of these 
inhomogeneities, 



2D(<j))5t g = V2DH-HIH 



(98) 



D. Conditions for self-reproduction 

Qualitatively, self-reproduction is present if the typi- 
cal change in the averaged field due to random quantum 
jumps ("diffusion") during one Hubble time is much larger 
than the deterministic change due to the classical equa- 
tion of motion, 



y / 2D(<p)5t > v((/))8t, St = H-\ 



(94) 



It is straightforward to see that this condition generically 
holds in the slow- roll regime. Using Eq. (|92() . we find 



(vStf 
DSt 



DH 



H 2 



as 4> — ► oo, 



(95) 



should not lead to a significant variation of H within a 
Hubble-sized region. (Otherwise, we cannot approximate 
the spacetime as a locally homogeneous de Sitter with a 
well-defined value of H.) Thus we obtain the condition 

^A 54>g = EdL^—^ §4> = ^DH^hicJ 1 ^ « 1, 
H H H 

(99) 

which is stronger than the condition Ij97(l by the factor 
\/ln cj 1 . However, this factor is an extremely slowly 
growing function of 0, which we may expect to contribute 
at most a factor of 2 or 3, and definitely no more than 
one order of magnitude. Therefore, we may ignore this 
factor and concentrate on the condition (|97ll . 
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To investigate whether the condition (|97|) may be vio- 
lated at large where i>*(</>) « «o, we can use the analytic 
approximations derived above. Using the formulas 1)00( 1 . 
(E2J), and we find 



V. ANALYSIS OF FOKKER-PLANCK 
EQUATIONS 

A. Self-adjoint form of the FP equation 



V2DH- 1 



4vrV Q'o 



3/4' 



AC 7 / 4 
4-7T 



(-Qo) 1/8 (Q^') 1/4 [3h|] 3/4 . 



(101) 



The constant coefficient b\ is small since it contains a 
large power of the inverse Planck mass, while other rel- 
evant energy scales are presumably much lower. Since 
(In if) , = (In if) the condition JHU) yields 



h (i 



— a/4 
KaVk) < i 



(102) 



Assuming the asymptotic behavior K(<fi) ~ </> s for large 
0, where s > 0, we simplify the above condition to 



6l</>8 S 4 < 1. 



(103) 



This condition is satisfied for large <f> if < s < | but 



leads to an upper bound, 



< 



s > 



In the latter case, for 6 > 



on values of <f) if 
v the fluctuations 



are so large that a single "quantum jump" (within a Hub- 
ble time 8t ~ PI 1 ) may bring the value of </> outside of 
the slow-roll range, e.g. reach the end of inflation. It is 
clear that the diffusion approximation breaks down for 
4> > <t>max- In other words, regions of the universe with 
4> > 0max cannot be described within the present semi- 
classical framework because quantum fluctuations are too 
large. A boundary condition needs to be imposed at 
4> — 0max if we wish to proceed using the diffusion ap- 
proach. One may impose an absorbing boundary condi- 
tion, arguing that regions with <f> > </> max disappear into 
a "sea of eternal randomness." Alternatively, one may 
impose a reflecting boundary condition, arguing that the 
"sea of randomness" emits and absorbs equally many re- 
gions. Since this boundary condition serves only to val- 
idate the diffusion approach, results can be trusted only 
if they are insensitive to the chosen type of the boundary 
condition. 

Models with K{(j>) ~ cj) s and < s < | do not exhibit 
the problem described above, however it is still necessary 
to impose a boundary condition at some value <p = </>max- 
For instance, the energy density for large <j> may reach 
the Planck value, which determines a (model-dependent) 
Planck boundary </> max such that K(4> max ) \Qq\ ~ Mp t . 
Results of applying the diffusion approach should be in- 
sensitive to the type of the boundary condition as well as 
to the precise value of <fi max . 



As is long known, the FP equation 1(11)) can be re- 
duced to a manifestly self-adjoint form, which is formally 
similar to a Schrodinger equation for a particle in a onc- 
dimensional potential. One can then use standard results 
about the existence of bound states, which correspond to 
stationary solutions of the FP equation. 

A stationary solution of the FP equation is of the form 



PpOM) = P{\)i<t>) 



(104) 



where P<\) (</>) is a stationary probability distribution. 
The function Pt\\ (</>) must be everywhere positive, in- 
tegrable, and satisfy the stationary FP equation, 



L P P W = (£>P ( a)) !00 - (v.P w ) + 3HP W 



XP 



(105) 

The proper time t may be replaced by another time vari- 
able, dr — T(<f))dt, where T((j)) is an arbitrary posi- 
tive function of cj>, and it is implied that the value of 
r is obtained by integrating J T(<p)dt along comoving 
worldlines [2i|. The FP equation in the time gauge r 
has the same form as Eq. ((105)) . except that the coeffi- 
cients D,v,H are divided by T(<f>). For instance, with 
T(4>) = H{4>) we obtain the "scale factor" or "e-folding" 
time variable r = j Hdt = In a, and the stationary FP 
equation for P^^(</>)is 



D 



,<M> 



v * p(sf) 



op(sf) _ , p(sf) 
(A) ~ Ai (A) ' 



(106) 



Self-reproduction is eternal if there exists a stationary 
solution Pf\) of Eq. I(105() with a positive A. This condi- 
tion is independent of the time gauge within the class of 
gauges dr = T((j>)dt [3. 

To investigate the existence of positive eigenvalues of 
the operator L p , it is convenient to transform the station- 
ary FP equation into an explicitly self-adjoint form. We 
perform the calculation following Ref. 23] and introduce 
new variables as follows, 



x = x(4>) = 



0,1, 



1 



P (A) (,/)) = ^(x)#- 3 / 4 ((/>) exp 



D 



d. 



v*(</>) 
D(<f>)' 



(107) 
(108) 



where the replacement <fi — > <f>(x) is implied in the last 
line. Note that P\{4>) is being multiplied by a strictly 
positive function of </>, so the positivity of ip(x) will entail 
the positivity of Pm j ((f)) . Under the replacements ((107(1 - 
1)108(1 . Eq. ((105(1 is transformed into 



ip,xic - U(x)tp = Xljj, 



(109) 
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where U(x) is the "potential" defined as a function of 
through the relation 



U{x) = -3H(<P) + 



3 (Drf) 



16 D 



2D 



v*,4, vj 

2 AD' 
(110) 

Equation (|109fl is the self-adjoint form of the stationary 
FP equation and is formally equivalent to a stationary 
Schrodinger equation for a particle in a one-dimensional 
"potential" U(x) with the "energy" E = —A. It is well- 
known that the ground state of a Schrodinger equation 
can be chosen as a nonnegative wavefunction. If we show 
that Eq. (|109J) has a ground state ip(x) with a negative 
value of "energy" Eq < 0, it will follow that there ex- 
ists a stationary probability distribution -P(A max )(</') cor- 
responding to the largest eigenvalue A max > of the FP 
equation. Then the existence of eternal self-reproduction 
will be established. 

In order to find out whether the ground state energy 
in a given potential U (x) is negative, it is convenient to 
use the variational principle, 

/ $ {-^,xx + Uip) dx 



Eq 



mm ■ 

i>{x) 



(111) 



J ifitpdx 

We may now substitute any test function ip{x) and ob- 
tain an upper bound on Eq. If we find such tp(x) that 
the integral in the numerator of Eq. (|lllf> is negative, it 
will follow that Eq < 0. It is sufficient to consider real- 
valued ip(x) that vanish at the boundaries of the range of 
integration. Integrating by parts and omitting boundary 
terms, we find 



il){-il), xx + Uil))dx= / (^) 2 + C/^ 2 dx. (112) 



If U(x) is everywhere positive, then clearly Eq > 0, so the 
only possibility to have a negative Eq is to have a range of 
x where U(x) < 0. Let us suppose that U(x) < within 
a range X\ < x < x 2 , with a typical value U(x) ~ — Uq < 
in the middle of the range. Then we may choose a 
test function tp(x) so that tp(x) — outside of the range 
[a;i,a;2] an d ip{%) ~ 1 within that range. The typical 



value of ip'ix) will be of order (x2 
we obtain a bound 

J (i// 2 + U4> 2 ) dx r 1 - UqI 



Xi 



I 1 , and so 



Eq < 



J ip 2 dx 



I 



-Uq- 



i. (113) 



The desired negative bound will be achieved if the width 
I of the interval is sufficiently large so that I 2 
However, the width I is constrained by 

Uq 



I < 



\U'(x)\> 



(114) 



Ynax x£[x 1 ,x 2 ] 

since the function U{x) must vary between and — Uq 
within the interval of width I. Hence, 

1 1 2 



Eq < -U 







\U'(x)\ 



u 







(115) 



In the next section, 
point x where U (x) - 



we shall use this condition near a 
= -Uq < 0. 



B. Stationary solutions in fc-inflation 

Let us now analyze a model of fc-inflation with the 
Lagrangian (|12|) . where 



K(4>) ~ 4> s , s> o, 



(116) 



for large </>. We assume that the conditions for the ex- 
istence of an inflationary attractor are met. Let us now 
estimate the "potential" U(x) at values of x that corre- 
spond to large <p. (Note that an infinite range of <j) may 
be mapped into a finite range of a;.) It is more conve- 
nient to analyze the behavior of the "potential" U as a 
function of (j). We find the following asymptotic behavior 
of parameters at <j) — ► oo: 

H(4>) - s/2 , v(4>) = const, (117) 
c,(0)~^-*-*, (118) 

(119) 



The dominant negative term in U{x) is —3H((f>) ~ —<j) s / 2 , 
while the dominant positive term could be only 



3 {D. 



D 



n (4 — n) 4> 



n-2 



16 D 4 
However, the condition (|97J) . rewritten as 

H, 



sfDH^-^- « 1, 



(120) 



(121) 



together with the power-law dependences of D((f>), H{<f), 
yields 





D H H 














H 


H (h 2 



H,4> 
IT 



< 1. 



(122) 
(123) 



Since the condition (|121f> must be satisfied for the entire 
range of 4> under consideration, the term l|120|) is negligi- 
ble compared with —3H. Similarly, the other terms de- 
cay as 0" 1 or faster for large <f>; for instance, v 2 /D <C H 
within the entire slow-roll regime (see Sec. lIVDll . There- 
fore, the potential U(x) can be estimated as U(x) w 
—ZH{4>) for most of the allowed range of (j>, excluding 
only a narrow range near the end of inflation. 

We now use the condition (|115fl to show that the 
ground state energy corresponding to the potential U(x) 
is negative. Since dx/d<j> > 0, and since the function 
H{4>) ~ s / 2 grows monotonically, the largest value of 
H x is at the largest allowed <fi. The bound l|115|) applied 
to an interval [4>\ , cj>2\ , upon using dx = \fDdtf), yields 

E < -3H(<t> 2 ) + D(<f> 2 ) 

w -3H(<h) + D{dp 2 )^ 2 « -ZH{cp 2 ) < 0, (124) 
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since D(j>~ 2 <§C H according to Eq. (|123[) . Thus we have 
shown that the ground state has negative "energy." This 
proves that the largest eigenvalue A max of the operator 
L p is positive, and hence that eternal self-reproduction is 
present. 

The bound Eq < —ZH{4>2), where 02 is essentially any 
value within the allowed range <p sr < 4> < </> max , means 
that the largest eigenvalue A max is approximately equal 
to 3P max , where PT max is the largest accessible value of 
H ((/)). It follows that the fractal dimension of the inflat- 
ing domain (defined in a gauge-invariant way in Ref . [T^ ) 
is generically close to 3. 

Finally, we note that the potential U(x) rj —3H(cf>) 
has a global minimum near the upper boundary </> max . 
It follows that the ground state ip(x) will have a sharp 
maximum near the upper boundary. Therefore, the sta- 
tionary distribution P(A max )(0) indicates that most of the 
3-volume contains <f> ~ m ax- However, it has been long 
known that that the distributions P c (</>, i) and P p (<j),t) 
are sensitive to the choice o f eq ual-time hypersurfaces, or 
the "time gauge" (see e.g. Efl Q Hi El E^). For 
instance, the 3-volume may be dominated by cf> ~ (f> max 
in one gauge but not in another; even the exponential 
growth of the 3-volume of the inflating domain is gauge- 
dependent 0. Nevertheless, the distributions P c ((j>,t) 
and Pp((j>,t) provide a useful qualitative picture of the 
global features of the spacetime during eternal inflation. 



C. Eigenvalues of the comoving FP equation 

We now consider the comoving distribution P c (</>, i), 
which describes the values of <f> along a single, randomly 
chosen comoving worldline. This distribution satisfies 
Eq. ©, 



dPc 
dt 



Lr.Pr. 



(125) 



which can be brought into a self-adjoint form similarly 
to the proper- volume FP equation. However, let us show 
directly that all the eigenvalues of the operator L c (with 
appropriate boundary conditions) are negative. It will 
follow that the late-time asymptotic of Eq. (|125|) is 



Pr 



,<) =^(A„)(0)e Aot , A o <0. (126) 

In other words, the probability P c (4>,t) tends uniformly 
to zero at late times, which is interpreted to mean that 
the evolution of <j> along any particular comoving world- 
line will almost surely (with probability 1) eventually ar- 
rive at one of the boundaries, either at c/>e or at </> max . 
Equation (|125fl can be written as a "conservation law," 



d t P c = d J, J{4>) = d^(DP c ) - u*P c 



(127) 



where J plays the role of the "current." The rele- 
vant types of boundary conditions, which were discussed 
above, are conveniently expressed in terms of the quan- 
tity J. Namely, the reflecting condition is J = 0, the 



absorbing condition is P c = 0, and the "exit-only" con- 
dition (to be imposed at the end-of-inflation point <Pe) 
is J = — u*P c , meaning that diffusion cannot bring the 
field from <f> — cj>E back into the inflationary range. To 
avoid considering the absorbing and the reflecting con- 
ditions separately, we shall impose at (f> = m ax a for- 
mal combination of an absorbing and reflecting condi- 
tion, J + aP c — 0, where a > 0. Suppose P(A)(0) is an 

eigenfunction satisfying L c P^x) 
ary conditions 



\P(\) with the bound- 



u 



(A) 



U))l 



= 0, 



J, 



(A) 



Pw)\ 

J (x) = d*(DP w )-v*P w . 



= 0, 
(128) 
(129) 



Note that the operator L c is self-adjoint with respect to 
the scalar product in the u ip n space, 



(Pi,P 2 ) = J Mx)Mx)dx= I Pi(0)P 2 (</))M(0), 
M(4>) = D((j>) exp 







-J 





(130) 

So let us consider the scalar product of Pr\) and L C P^, 

I[P W ] = (P (X) ,L C P W ) = y 0maX P( A ) (LcP(x)) Mdct>. 

(131) 

By assumption, L c P(a) = ^P(X) , thus 



I[P(X)] = x 



Pf x) Md(j>. 



(132) 



Since M(cf>) > 0, we will prove that A < if we show that 
I[P(X)\ < ; 

Writing L c P(a) = d$ J(a) and integrating Eq. I|131|) by 
parts, we find 



I[P W ] = MP W J W \% 



#J (A) ^(MP (A) ). 



(133) 

Using the boundary conditions l|128|) . we can estimate 
the boundary term MP^x)J(X) as follows, 



(^W(A))l;r 

= -aMP ( 2 A) (0 max 



) + v*MP 2 {x) {<j )E ) <0 (134) 

because v*(0) < and a > 0. It remains to put a 
bound on the integral in Eq. (|133fl . It is easy to see 
that <9^(M P(a)) is proportional to J(\), namely 



(135) 



Since M(<f>)/D(<f>) > 0, the integral in Eq. I|133|l must be 
negative (note that J^x){<p) is not everywhere zero): 



dcf>J w d^MP w ) 



, M 



(136) 
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Thus we have shown that -f[-P(A)] < 0, which proves that 
every eigenvalue A of L c is negative. 

D. Exit probability in fc-inflation 

We have seen that the evolution of the averaged field 
(f> is heavily influenced by random "jumps," which can 
even bring the value of <j> to the upper boundary <fi max 
where the semiclassical description breaks down. Thus, 
we may trace the evolution of (j> along a single comov- 
ing trajectory, starting from an initial value 4>a, and ask 
for the probability p C xit of eventually exiting through the 
end-of-inflation boundary (f>E while staying away from 
the upper boundary </> m ax- In this section, we show that 
this exit probability is approximately equal to 1 if 4>q is 
sufficiently far away from </> max - 

The exit probability can be found as follows (a similar 
method was used in Ref. H^)- Let us assume that the 
initial distribution is concentrated at = 00: i-e. 



Pr 



,t = Q) = 6(<l>-<M, 



(137) 



where 4>e < O < 0max> and that the distribution P c ((j), t) 
is known at all times. Due to the conservation law (|127|l . 
the probability of exiting inflation through = (f>E during 
a time interval [t, t + dt] is 



dpexit = Jdt = {d<j, (DP C ) - v*P c 

= -V*{4>E)Pc(4>E,t)dt 



dt 



(138) 



(note that u*(</>_e) < 0), hence the total probability of 
exiting through <p = 4>e at any time is 



rfpexit = -v*{4>e) I P c (4>E,t)dt. (139) 



Poxit 



To compute p ex it directly without knowledge of P c (<f>,t), 
we define an auxiliary function 



Q{<t>) 



diP c (<M), 



(140) 



so that 

Poxit = -v*(<Pe)Q(<Pe)- (141) 
It is easy to see that the function Q(4>) is a solution of 

L C Q{4>) = P c (<t>, t = oo) - P (<j>, t = 0) = -5{4> - fa), 

(142) 

since P c (4>, t — oo) = 0. (Note that the function 
v*((f>)Q{4>) satisfies a gauge-invariant equation, 



84, \—{v*Q) ) -v*Q 



-8{4>-4>o), (143) 



which reflects the fact that p C xit is a gauge-invariant 
quantity.) The boundary conditions are the time- 
integrated conditions 11281) . 



8<t> {DQ)\, = 0, Q(0 max ) = O, 



(144) 



where we have chosen the purely absorbing boundary 
condition at (j> = ^ max because we are now interested 
in regions that never reach that boundary. Rewriting 
Eq. (Jig) as 

3$ [dj, (DQ) - v*Q] = -5{<t> - <p ), (145) 

we may immediately integrate, 

84, (DQ) - v*Q = C t - 6(<j> - <f> ), (146) 

where we note that G\ = p C xit due to the boundary con- 
dition at <fi — <pE- Then we obtain the general solution 



Q(<f>) = c 2 



1 



D(<f>) 



exp 



D(4>) 



i/(Ci-e(4/-<M)exp 



, d' 



(147) 



The constants of integration Ci,2 are determined from 
the boundary conditions ljl44|l . and the final result is 



Poxit 



R(4>E 

R((j>) = exp 



r<t> 


max y 


Jd> 





(148) 
(149) 



The above equations are valid for any i>*(</>) and D(<j>); 
we shall now specialize to the case of fc-inflation. Since 
v* < and is approximately constant, while D(<p) > 
and grows with 4>, the function R((f>) is approximately 
equal to 1 within a certain (model-dependent) range of 
0, namely for tfin < <f> < </> ma x, where (j>R is determined 
by the condition 



1 



D' 



DUr 



(150) 



Since R(<fi) quickly approaches zero for <f> < <f>R, we may 
neglect R(4>e) "C 1 and express the exit probability as 



Pexit — C\ 



st- Rim' 



It follows that 

Pcxit 
Pexit 



»o < <pr\ 
- 4>o 



OR 



»0 > 4>R- 



(151) 

(152) 
(153) 



A worldline starting in the middle of the range [4>r, <^ ma x] 
will exit either at <f> — <f>E or at the upper boundary 
— ^max with nearly equal probability. Therefore, we 
may interpret the range [4>r, <fimax] as the "runaway dif- 
fusion" regime. A typical comoving worldline will surely 
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(with probability p cx it ~ 1) reach the reheating boundary 
(j)E only if the initial value 0o is outside of the "runaway 
diffusion" regime, i.e. 0o < <Ar- The value 0_r is highly 
model-dependent, but generically 0# <C 0max- There- 
fore, the existence of the "runaway diffusion" regime lim- 
its the independence of fc-inflation models on initial con- 
ditions. 

Note that the "runaway diffusion" regime is absent in 
models of potential-driven inflation, where we have 

V H z 

, DU) = — , (154) 

3H' KV> 8tt 2 ' v ' 



ff(0) « kVV, u*(0) 



and the meaningful range of is < i^ max , where ifi, l;l .x 
is the Planck boundary determined by V(<p max ) 
One finds 

8^ / 1 1_ 

'3k 4 {v(cf)) ~ V(0 maj 



R(4>) = exp 



(155) 



so R(<f>) is negligibly small almost all the way until the 
boundary = (j) max . It follows that 0# ~ 0max, so the 
exit probability in potential-driven inflation is p C xit ~ 1 
for any 0o within the allowed range 0o < ma x- 
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Appendix A: ATTRACTOR SOLUTIONS IN 
POTENTIAL-DRIVEN INFLATION 



In this appendix, we derive the attractor solution in 
models of potential-driven inflation with a power-law po- 
tential V(4>) oc 0™. This calculation serves as an illustra- 
tion of the general method for analyzing attractors that 
we developed in Sec. IIIII While the existence of attrac- 
tor behavior in potential-driven inflation is long known, 
it has not been stressed that a unique attractor solution 
can be singled out in the phase plane (see e.g. Ref. 
where a statement is made to the contrary). 

We consider a cosmological model of inflation driven 
by a minimally coupled scalar field with a potential 
V (0). The equations of motion are 



+ 3ff0 + V'(0) = 0, 



U 2 + V(cp) 



3M 



pi 



8tt 



H 2 = —H 2 . 

K 



(Al) 
(A2) 



Considering ^ as a function of 0, i.e. = u(0), we can 
reduce the above equations to 



du{4>) _ V'{cf) 



?>K\I- + V{4>) = g{u,4>). (A3) 



For the purposes of further analysis, we assume that the 
potential V(<p) has a monotonically growing behavior at 
— > oo, going as V(<p) oc <fi n with n > 2. 



1. Determining the attractor solution 



According to the criterion developed in Sec. IIII Al at- 
tractors are determined by their behavior at — > oo. In 
the phase plane (0,0), an attractor solution = u*(0) 
should have the property that neighbor solutions grow 
significantly faster than u»(0) at — > oo. In the case 
of potential-driven inflation with growing V(0), we have 
(at fixed u) 



g(u, 0) ~ max(V',VV) 



oo as ' 



(A4) 



Hence, every solution u(0) grows at — > oo. How- 
ever, a generic trajectory grows exponentially fast for 
large 0, while we expect the attractor trajectory to grow 
slower than exponentially (e.g. polynomially). To make 
the attractor behavior apparent, we change variables as 
it(0) = F(0)u(0), where F(0) is a fixed function that 
will be determined below. The function F(<ft) should have 
polynomial growth, F(0) ~ s , such that the attractor 
solution is u(0) — * const at — > oo. 

After the change of variables, the equation for m(0) is 



du g(Fu, 0) - F'u 



F 



V 
F 2 u 



3k 



V 

J2 



EL 

~F 



u = g(u,(f>). (A5) 



A solution u(0) — > const at — * oo will exist if g{u 1 < 
has an "asymptotic root" u < such that 



lim g(u ,<j>) 



0. 



(A6) 



Since F'/F ~ 0" 1 , the last term in Eq. (|A5|) is always 
dominated by the second term at fixed u as — » oo. 
Therefore, an asymptotic root will exist if the first two 
terms cancel each other as — > oo. This requires that 
one of the two sets of conditions hold in the asymptotic 
limit — > oo : 



YL y y 

F 2 ~ F ' 



p2 



V 



» 1, 



V 



< 1. 



(A7) 



(A8) 



For a power-law potential, V(0) = A0", and F(<fr) — 
0™, it is straightforward to see that only the first set of 
conditions can be met, which yields 

n V 
m=--l, F(0)=0^- x ~— . (A9) 



The "asymptotic root" uq is 



u 



3k 



(AlO) 



(The value of Uq is negative since < for large 0.) With 
the choice F{<j>) = 0"/ 2_1 , the attractor is determined as 
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the (unique) solution ((/>) which approaches a constant 
as (j> — ► oo. Other solutions grow exponentially with <j). 
In the original variables, the attractor solution is 

u*{4>) ps F(<p)u ~ c^ -1 as ^ — * co. (All) 

This is the well-known slow-roll attractor behavior in 
potential-driven inflation, ~ V /W. 



This is the familiar slow-roll approximation. Higher- 
order asymptotic approximations may be determined by 
iterating Eq. (|A5(I . considered as an equation for u((f)) 
with a given du/d<f>. After some algebra, the next-order 
approximation to the attractor solution is found as 



2. Approximate expressions for the attractor 

A first approximation to the attractor solution is 
«*(</>) ps Uq = const, which corresponds (in the original 
variables) to the solution of the equation 



= ft = 



V 



3kVV 



(A12) 



V 



3kVV 



AW" - 3V 7 - 
36k 2 V 2 



(A13) 



As expected, this expression reproduces the slow-roll ex- 
pansion of Ref. [l7j up to terms of first order. 
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